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Abstract 



We have studied sidebranching induced by fluctuations in dendritic 
growth. The amplitude of sidebranching induced by internal (equilibrium) 
concentration fluctuations in the case of solidification with solutal diffusion is 
computed. This amplitude turns out to be significantly smaller than values 
reported in previous experiments. The effects of other possible sources of fluc- 
tuations (of an external origin) are examined by introducing non-conserved 
noise in a phase-field model. This reproduces the characteristics of sidebranch- 
ing found in experiments. Results also show that sidebranching induced by 
external noise is qualitatively similar to that of internal noise, and it is only 
distinguished by its amplitude. 

PACS numbers: 68.70,+w, 81.10.Aj, 64.70.Dv, 05.40.-a, 81.30.Fb, 05.70.Ln 
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I. INTRODUCTION 



Dendritic growth in nonequilibrium systems has been extensively studied during the last 
few years [|l|-|7|]. A feature which remains a main point of interest is the study of sidebranch- 
ing, which is the secondary branches that appear at both sides of the main dendrite. The 
question of how its frequency and amplitude are determined has not yet been fully solved. 
Two scenarios have been proposed to explain the origin of sidebranching. One of them 
states that periodic deterministic oscillations at the tip P-flCfl generate correlated branches 
on both flanks of the dendrite |||§. A possible source of these tip oscillations was suggested 



in Ref . [fTTJ , where it was argued that the spreading rate of the wavepacket that characterizes 
sidebranching might become large enough so that the tip could undergo oscillations or other 
instabilities. This is predicted to occur in the limit of small surface tension anisotropy. The 
other scenario proposes that sidebranching is due to selective amplification of fluctuations 
near the tip of the dendrite [|T2|-f2"7[| . In this case, branches appear to be mostly uncorrelated. 



In this paper we will study this second scenario by means of a phase-field model |]25|-f40 
and, in particular, we will focus on the issue of an external vs. internal origin of the noise. 

In a frame of reference moving with the tip of the dendrite, sidebranching can be seen 
as a wave that propagates along the dendrite away from the tip at the same velocity as 
the tip. An appropriate characterization is provided by the amplitude and wavelength. 
Barber et al. [[14] studied the evolution of time- dependent deformations of the needle crystal 
(Ivantsov) solution of the two-dimensional symmetric model of solidification in the limit 
of small Peclet number within a WKB approximation. The amplitude of a localized wave 
packet grows exponentially as z 1//4 , where z is measured from the tip along the symmetry 
axis of the dendrite, as the packet moves down, provided the initial packet contains modes 
of arbitrarily small frequencies. Moreover, the wavelength of the packet increases as z 1//4 . 
Pieters jn| obtained the same amplitude and wavelength dependence on z as in Ref. [13 



both analytically and by numerical integration of the boundary-layer model. Langer f!6 
concluded, from a similar analysis to that of Ref. |14j but performed in three dimensions, 
that noise of some kind can be the origin of sidebranching, but that thermal fluctuations are 
not strong enough to entirely explain the phenomena. However, more recently, Brener and 
Temkin |L7|| considered anisotropic needle crystals in three dimensions and concluded that 
experimentally observed sidebranching could be explained by considering noise of a thermal 
origin. The growth of the sidebranching amplitude was found to behave exponentially as a 
function of z 2//5 , which is faster than the z 1 ^ dependence obtained in the axisymmetric case 
T6|j . The sidebranching wavelength was found to be a function of z 1//5 , very similar to that 



obtained in the axisymmetric case. Dougherty et al. [18[] studied sidebranching in NH 4 Br 
dendrites, where rather uncorrelated variations in phase and amplitude of the branches were 
observed. They determined the amplitude of the sidebranching and its mean wavelength 
by looking at the power spectrum of the data obtained by measuring the half-width of 
the dendrite at a fixed distance z from the tip at different times. The behaviour of the 
amplitude was qualitatively like that predicted in Ref. | up to a certain value of z, after 
which the linear theory is presumably no longer valid. An equivalent exponential growth of 
the amplitude with s 1 / 4 , where s is the arclength, was also found in Ref. However, in 
Ref. [|18j no variation of the mean frequency in the spectral peak was obtained for different 
z. Finally, Dougherty et al. also observed that side branches separated by more than about 
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six times the mean wavelength were uncorrelated. Weak correlation between opposite sides 



of the dendrite when no external forcing was applied to it was also found in Ref. |20|| . The 
common feature of all these experiments is that sidebranching appears to be due to the 
selective amplification of natural noise and not to the existence of some intrinsic oscillation 
or limit cycle. 

Bisang and Bilgram 1 211,221 found quantitative agreement between the predictions for the 
linear regime in W7\ and their results in experiments on xenon dendrites in three dimensions. 
They concluded that Brener and Temkin's theory correctly describes the sidebranching be- 
haviour of dendrites for any pure substance with cubic symmetry and thus thermal noise 
was concluded to be the origin of the sidebranching observed in their experiments. 

In the last decade there has been an increasing use of phase-field models to deal with 
dendritic growth problems. They are a very useful and practical tool to simulate such kind 
of processes and a good alternative to the integrodifferential equation which can be derived 
from the classical sharp-interface model. In the phase-field models an order parameter 
or phase field is defined, which avoids tracking the interface and naturally includes the 
physical boundary conditions at the interface. 

Up to now, few studies of sidebranching phenomena with the phase-field model have been 
carried out. It has been shown that the inclusion of a noise source in the phase-field model 
equations enhances the emergence of uncorrelated sidebranching |28[] without affecting the 
velocity and radius of the tip |p9[ . Moreover, when the dendrite tip is periodically forced, 



the sidebranching appears to be correlated at both sides of the dendrite |33| as has been 
observed in some experiments p0| , |2^| . In particular, sidebranching can be regulated by 
spatially homogeneous time-periodic variation of the melting point induced by oscillations 
in the external pressure or by periodic heating generated by a dissipative electric current 

The deepest insight into the study of sidebranching with a phase-field model have been 



carried out recently in Ref. ||26|| . They included thermal noise in a two-dimensional phase- 
field model of solidification controlled by heat-diffusion in a way which was consistent with 
both bulk and interfacial equilibrium fluctuations, as has been done previously with the 
sharp- interface model equations P^J13[| . Karma and Rappel |26[ obtained good quantitative 
agreement between the computed sidebranching amplitude as a function of distance to the 
tip and the prediction of the linear WKB theory for anisotropic crystals in two dimensions. 
Sidebranching wavelength very close to the tip was found to increase with z faster than pre- 
dicted by the WKB theory, but this could be explained after considering that perturbations 
generally get stretched as they travel along the sides of curved fronts. Further from the tip, 
the value of the sidebranching wavelength saturates. 

Although there is general agreement in that thermal (internal) fluctuations are enough 
to explain the amplitude of the dendritic sidebranching, one should be aware that evidence 
along these lines has been achieved in experiments of heat-controlled solidification of pure 
substances. As up to now there is a lack of predictions of sidebranching amplitudes for 
solutal dendrites, experiments of these dendrites can only show qualitative agreement with 
theoretical results. In this paper we address the question of sidebranching characteristics in 
the presence of external vs. internal fluctuations. First of all, we obtain a prediction of the 
effects of internal noise on sidebranching amplitudes for solutal dendrites. A comparison of 



the theory with available quantitative experimental results 1T3] shows that there are seri 
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ous indications that in some experiments internal thermodynamical fluctuations could not 
account for observed sidebranching activity. In this case some other source of fluctuations, 
of an external origin, should be called on. Some of the consequences derived from adding 
a non-conserved noise source in a two-dimensional phase-field model are examined. This 
noise is of a different nature to what one should employ to provide an account of internal 
fluctuations |25|- ^7|^2| , ^3| . However, our simulations qualitatively reproduce the properties 
of the noise-induced sidebranching derived analytically [0,^] and observed experimentally 
18| , l9[] . We conclude that although thermal noise is not always the main origin of dendritic 



sidebranching, its qualitative characteristics are common for noise-induced sidebranching 
independently of its origin. A detailed quantitative study of sidebranching activity could 
therefore be useful to elucidate the origin of the noise in specific experiments. 

In Section II we predict the sidebranching amplitudes for solutal dendrites with thermal 
fluctuations. In Section III the model and the numerical method used to solve the equations 
are described. Numerical results on the behaviour of the sidebranching characteristics as 
well as comparison with theoretical predictions and experimental results are described in 
Section IV. Conclusions derived from these results are outlined in Section V. 



II. SIDEBRANCHING AMPLITUDE IN SOLUTAL SOLIDIFICATION 

Available theoretical predictions on sidebranching amplitudes have been formulated for 
dendrites grown from a pure substance and controlled by heat diffusion. Here we will 
consider dendrites appearing in isothermal growth of mixtures controlled by diffusion of the 
solute. We start from the Langevin formalism for solidification due to Karma |4^j4*3[1 . In this 



formalism the usual sharp-interface model for solidification is completed with noise terms 
constructed with the requirement that they give the correct bulk and interfacial equilibrium 
fluctuations. The resulting diffusion equation for a mixture in isothermal conditions is: 

BO 

-^ = W-V-#,t), (2.1) 

with the following boundary conditions at the interface: 

(C L - C s )v n = n • [D S VC S - D L VC L ] + n ■ [q L - q s ] (2.2) 

m E C L + Tk + ^=T m -T + V (r, t), (2.3) 
H 

where v = S,L denotes the phase, C u is the concentration, Tm — T is the undercooling, 
D u is the diffusion coefficient, m E is the absolute value of the (negative) T{Cl) slope of 
the coexistence curve, and k, n and v n are the curvature, the normal unitary vector and 
the normal velocity of the interface, respectively. V = aT/L, a is the surface energy and L 
the latent heat per unit volume, q and 77 are random forces whose statistical properties are 
given by 

ti(r,tW(r',t')) = 2D„C„(r,t)S(r,r')S(t-t% j (2.4) 
1 ./l + |V_ L £(r L ,t)| 
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where the interface is parametrized by the vector r_i_ + £n. 

These equations can be mapped under several approximations into the corresponding 
Langevin model for free solidification of a pure substance 0,^,^]. First we assume a 
constant concentration gap in the mass conservation Eq. |2.2| , i.e. Cl — Cs = AC pa AC eg , 
the value corresponding to equilibrium at temperature T. This is in principle valid for 
small curvatures and velocities. A similar approximation is assumed in the intensity of the 
bulk noise q%(r,t), substituting C„(r, t) by the equilibrium value C* q in Eq. |2.4j . Further- 
more in the intensity of the interfacial noise rj^(r±,t) we employ the Clausius-Clapeyron 



equation for dilute alloys |43| to make the substitution ~ ~A<pf m Eq- Within 
these approximations the process of isothermal solidification of an alloy is equivalent (in- 
cluding thermodynamical fluctuations) to the (heat diffusion controlled) solidification of a 
pure substance, whose specific heat, latent heat and melting temperature are given by 



k 



B 



(AC**)' 

icq 



(C e ^ 3 



T 
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AC e i 
(Cl q f 
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m E 



where the diffusion field is now a temperature field given by 

T(T,t)-T M =(C v -C?)m E 



(2.6) 

(2.7) 
(2.8) 

(2.9) 



and the rest of parameters remain unchanged. This can be checked by direct substitution 
in the Langevin equations. Therefore, the sidebranching induced by thermodynamical fluc- 
tuations should be the same in both situations. Sidebranching amplitude is predicted to 



depend on the distance z along the dendrite clXIS clS 

H7T 



A(z) = pS'exp 



4( z ) 

3cr*zp 2 



1/2N 



(2.10) 



where p is the tip radius, x (z) is the shape of the dendrite, and the operating mode of the 
dendrite is given by the parameter a* defined by 

a* = 2D d /p 2 V } (2.11) 

where V is the selected velocity and do is the capillary length. The dimensionless noise 



amplitude, for a (i-dimensional thermal dendrite, is then known to be [17.26 



£3 _ 2k B T 2 M c v D 



(2.12) 



L 2 p l + d V ' 

Applying the mapping above, the corresponding result for a d- dimensional solutal dendrite 



is 



2C e L q D 



(AC e if p 1+ "" 



V 



(2.13) 
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This result will be used below to compare the prediction for the case of internal noise with 
experiments on solutal dendrites, for which there are not many quantitative experimental 
results available. We focus on the experiments performed in Ref. []TH| with ammonium 
bromide dendrites growing from supersaturated aqueous solution in isothermal conditions. 
In this experiment the precipitate front advances by incorporating solute particles instead of 
rejecting them, which makes it differ by several details from standard solutal solidification. 
Nevertheless, the above result (Eq. |2.13| ) can be obtained by slightly adapting the performed 
mapping. In this case the system is on the high concentration side of the phase diagram of 
the mixture, for which Tm is that of the solvent, no longer close to the temperature of the 
experiment. In this case it is convenient to write the Gibbs-Thomson equation as 



mE{Coo -C l ) + Tk+^ = T s -T + V {r, t), 



(2.14) 



where is the concentration of the dilution, T$ is the saturation temperature for that 



concentration, and tue is now the (positive) T{C£) slope. The same results of Eqs. |2.10| and 



2.13 are obtained by applying the mapping of Eqs. |2.6|, 2.7 and 2.8 with a diffusion field 



T(r,t)-T M = (CZ«-C u )m E , 



(2.15) 



where we have used the relation Ts — T = ra^Cco — C e ^). 

Now we compare the prediction of Eqs. |2.10| and 2.13 with the experimental results on 
supersaturated solutions of Ref. ||18|| . These experiments were performed at a supersatu- 
ration A = (C^ — C eq )/{Cs — C eq ) = 0.007 and a saturation temperature of 56°C. The 
characteristics of the selected dendrite are p = 4.0/xm, V = lAifim/s and a* = 0.081. We 
have employed a value of D = 2.6 x 10 _5 cm 2 /s and the values corresponding to a tem- 
perature of 100°C for the equilibrium concentrations: C e ^ = 0.99 x l^Nxmolec/m 2, and 
Cg q = 2.48 x 10 4 iV \molecj 'm? . Since this temperature is much higher than that of the exper- 
iment, the resulting value S = 5.96 x 10~ 5 constitutes an overestimation of the theoretical 
value. 

In Fig. 1 we show the calculated amplitude A(z) of the sidebranching induced by internal 
fluctuations for this value of S. We consider the theoretically predicted shape x(z) = (^z) 3 ^ 5 , 



as was considered in Ref. |L7[]. The result corresponding to the actual temperature of the 
experiment would be placed below the represented curve. In the same figure we plot the 
experimental results of Ref. [pl|. We see that experimental amplitudes are approximately 
one order of magnitude larger than the overestimated theoretical values. 

Therefore the predicted amplitude of the sidebranching when it is due to statistical 
noise is, for the experiments of Ref. [[18] 
amplitude experimentally observed. Thus, 



at least one order of magnitude smaller than the 
we are led to conclude that thermodynamical 
fluctuations are not enough to explain the sidebranching amplitude in some experiments. In 
this estimation we have assumed three-dimensional dendrites even though the experiments 
of Ref. |18| are intended to be quasi two-dimensional. The analogous calculations in two 
dimensions show an even lower sidebranching amplitude, i.e. a larger discrepancy with the 
experiments. 

There is a shortcoming in the predictions above when applied to supersaturated ex- 
periments. By their own nature, supersaturated dilutions are not in the diluted limit, as 
assumed in the theoretical analysis ||43|| . In the experiments of Ref. ]18[ the concentration 
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is as high as 16% of solute molecules. Indeed, the whole Langevin formalism is constructed 
in order to guarantee that the concentration fluctuations in a small volume AV is 



(Aa) 2 ) = (2.16) 
which is the equilibrium value for diluted solutions. In fact, for a concentrated solution Eq. 



2.16| should be replaced by 



\dcJ PtT 

Since the sidebranching amplitude in supersaturated solutions has been found to be at least 
1 order of magnitude larger than predicted in the diluted approximation, one concludes that 
it would be necessary that the derivative of the chemical potential is 2 orders of magnitude 
greater than that given by the diluted approximation |^ ~ ^ in order to explain the 
experiments by internal noise. Therefore, most likely this internal noise is really not strong 
enough to account for the observed sidebranching. As we are not aware of quantitative 
data on thermodynamical properties of supersaturated solutions that would permit us to 
improve the estimations above, a definitive answer on the amplitude of sidebranching in 
these dendrites remains open. In any case these results call for experimental quantitative 
measurements in solutal dendrites grown in diluted conditions, where Eq. 2.13| properly 
applies. 



III. MODEL AND NUMERICAL PROCEDURE 

We have performed simulations of dendritic growth by employing a phase-field model 
for solidification. In this model both phases and their interface are treated indistinctly, and 
discriminated by an effective non-conserved order parameter or phase field <f), which takes 
different values in each phase (0 and 1 in our simulations). This field changes smoothly across 
an interface region of finite thickness, and its dynamics is coupled to that of the diffusion 
field in such a way that the sharp-interface model is recovered in the limit of vanishing 
interface thickness, controlled by a new small parameter e. The equations of the model read 
explicitly : 



— e 



d_ 

dx 



e 2 r{9) 



Of 



dy 



= 0(1- 

Ve 2 — 
dy 



1 



+ 3Oe/?Aw0(l 



+ e 2 V rf{9)V(t) 



(3.1) 



6O0 3 + 3O0 4 ) = V 2 u + ip(x, y, t) 
' at 



(3.2) 



where u(r, t) is the diffusion field and A is the dimensionless undercooling. Lengths are scaled 
in some arbitrary reference length lu, while times are scaled by uj 2 /D. In these equations 9 
is the angle between the x-axis and the gradient of the phase field. rj(9) is the anisotropy 
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of the surface tension. The anisotropy of the kinetic term is then given by t(9)/t](9). (3 is 
equal to y^p and d is the capillary length. 

The external noise is introduced through the additive term ip in the heat equation. This 
choice is not unique and is justified here only for simplicity. Because of its external origin, 
the noise is not assumed to satisfy a fluctuation-dissipation relation. Furthermore, it is 
generally assumed to be non-conserved, as opposed to the more usual case of thermal noise, 
which would enter the model equations as a stochastic current (i.e. conserved noise) in the 
heat equation, and an additional stochastic term in the phase- field equation f2~6| . In our 
simulations the noise term is evaluated at each uncorrelated cell of lateral size Ax simply as 
I ■ r, where / denotes the amplitude of the noise, and r is a uniform random number in the 
interval [—0.5, 0.5]. The phase-field model equations have been solved on rectangular lattices 
using first-order finite differences on a uniform grid with mesh spacing Ax. An explicit time- 
differencing scheme has been used to solve the equation for 0, whereas for the u equation the 



alternating-direction implicit (ADI) method was chosen, which is unconditionally stable |29 



The kinetic term has been taken as isotropic, which leads to t{6) = mrj(8) with constant 
m. A four-fold free energy anisotropy a = o"(0)(l + 7cos(4#)) has been considered. This 
gives rise to dendrites growing with perpendicular side branches. The values of 7 taken 
were always smaller than 0.0625, which ensured that we obtained rounded shapes such as a 
parabola because forbidden directions were avoided |5|]. 

The growth morphologies were obtained by setting a small vertical seed (0 = 0, u = 0) in 
the centre of either of the two shortest sides of the system and imposing = 1 and u — — 1 
on the rest of the system. We have considered symmetric boundary conditions for and u 
on the four sides of the system, although they do not influence the results presented in this 
paper. 

We have used a set of phase-field model parameters that give rise to a growing needle 
without sidebranching for every anisotropy 7 considered when no noise is added to the 
simulations. This ensures that the sidebranching observed when noise is present is not due 
to computational rounding. The fixed phase-field model parameters for all the simulations 
are j3 = 400, m = 20 and e = 0.003. The values of A and 7 have been varied in the range 
0.5 — 0.7 and 0.045 — 0.060, respectively. The noise amplitude is varied in the range 5 — 150. 
The time and spatial discretizations were kept constant in all the simulations with values 
Ax = 0.01 and At = 10~ 4 . Two system sizes were used in the simulations. A system 
with 1000 x 1500 grid points was used to observe fully developed sidebranching including 
tertiary arms and to have enough statistics in order to compare with measurements of the 
sidebranching correlation at both sides of the dendrite. Additionally a 500 x 500 grid points 
system was used when the data sets did not require very extensive statistics. In Fig. 2 a 
growing dendrite is shown with 7 = 0.045, A = 0.6, I — 11 at a time t = 1.5. The velocity 
and the radius of the tip are very weakly affected when noise is introduced. However, side 
branches appear at both sides of the main dendrite, yielding approximately a 90° angle, 



as was observed in ||18|| . Further down the tip one can clearly observe competition between 
branches which gives rise to a coarsening effect. When branches reach the vertical boundaries 
of the system, they are stopped by the effect of the symmetric boundary conditions. This did 
not affect the measurements presented in this paper, where we have focused on the region 
between the tip and a point approximately 150 grid points down from the tip (grid points are 
marked on the axes of Fig. 2). This region corresponds approximately to that considered 



S 



in the experimental work in Ref. |18] and, moreover, it is an appropriate region in order 
to compare the behaviour of sidebranching with the analytical results obtained from linear 
perturbation theory. In the longer runs, and in order to avoid working with unnecessarily 
large systems, we have performed periodic shifts of the complete system. We have checked 
that this did not affect the results of the simulation. 



IV. SIDEBRANCHING CHARACTERISTICS 

In order to study the sidebranching induced by noise we have measured the half-width 
h z (t) of the dendrite at various distances z behind the tip as a function of time. In order to 
have a comparable amount of data as in Ref. |18| (they recorded around 240 oscillations of 
the amplitude for each z) we needed to simulate a dendrite four times longer (t = 6) than 
the one shown in Fig. 2. The half width of the dendrite as a function of time at a distance 
z = 40 grid points from the tip as well as its power spectrum P z (f) are shown in Fig. 3. 
The same type of data corresponding to a point further from the tip (z = 100 grid points) 
is shown in Fig. 4. The data used to compute the power spectra were six times the lengths 
shown. Both sets of data have a strong resemblance to those obtained experimentally |18 . 
We also computed the cross-correlation function 

C(t) =< [h Lz (t + 1) - h Lz ][h Rz (t) - h Rz ] > /a L a R , (4.1) 

where hi z (t), h Rz (t), <jl and a R are the half-width functions and their standard deviations 
for the two sides of the dendrite at the same distance from the tip. We found that C(0) is 
approximately 0.4 for points very close to the tip and that its value decreases very quickly 
to when increasing z. Moreover, in the points closer to the tip, C(t ) drops to zero after 
six oscillations, that is, the time to nucleate six side branches. The same behaviour was 



observed in the experiments in fL8fl . In simulations where we used a smaller data set, the 
values of C(t ) were larger than those found in the experiments, as was to be expected. 
Besides the predictions for the amplitude commented on above, the behaviour of the 



wavelength is another main feature of sidebranching. It depends very weakly on z fll7| , |26 
A ~ z 1 / 5 . Despite the fact that experiments p~8| , |T9l showed the predicted dependence of the 
amplitude on z, the same frequency of the spectral peak for different z was found. 

In order to observe the dependence of the sidebranching amplitude on z in the simu- 
lations, the phase-field model simulations were run on a 500 x 500 grid points system. In 
Fig. 5 we show the square root of the area under the spectral peak as a function of z. 
This representation gives the amplitude behaviour of the sidebranching as a function of the 
distance to the tip. The data were obtained with 7 = 0.045, A = 0.6, I = 19 until a time 
t = 0.5. The amplitude is found to increase exponentially with z 2 ^ for distances z < 100. 
Thus the behaviour of the data obtained in the simulations is consistent with the linear 
analysis p!7|j26|] up to a certain value of z. In []18| similar behaviour was found in the linear 



regime, although they found a saturation of the amplitude further down from the tip. This 
could be due to some boundary effects which stop the growth of the side branches. 

The dependence of the sidebranching wavelength on the distance to the tip was studied 
with the data recorded in the same simulation as for Fig. 5 and after performing the power 
spectrum. The wavelength was found to remain constant in the interval of considered z, 



which is consistent with the observations in f26 |. 



9 



V. CONCLUDING REMARKS 



We have studied the sidebranching induced by fluctuations in dendritic growth, for which 
there is good quantitative agreement between theory and experiment for thermal dendrites 
(i.e. controlled by heat diffusion). This agreement has only been qualitative for solutal den- 
drites, as there was no prediction available for their sidebranching amplitudes. In order to 
be able to perform a quantitative comparison, we have obtained an estimation of sidebranch- 
ing amplitudes originated by internal fluctuations in solutal diluted dendrites. The resulting 
values appear as much smaller (at least one order of magnitude) than those observed in some 
experiments (performed in concentrated solutions). This can be attributed to the effect of 
other noise sources, of an external origin. To confirm this conclusion, it would be necessary 
to make use of quantitative experimental results obtained in more diluted conditions. 

We have obtained the effects of non-internal fluctuations on dendritic sidebranching by 
introducing non-conserved noise in a phase-field model for solidification. Our simulations 
qualitatively reproduce well previous experimental and analytical results. In particular, we 
have reproduced the dependence of the sidebranching amplitude on the distance from the 
tip, confirming the validity of what was previously obtained with a linear theory including 
internal noise. Thus, the qualitative behaviour of the sidebranching amplitude, when this 
is due to the selective amplification of fluctuations, is basically independent of the origin 
of the noise. In conclusion, qualitative concordance between experimental results does not 
directly imply a common source of fluctuations and therefore, a careful quantitative study 
of sidebranching activity may help to elucidate the origin of the dominant noise in each 
experiment. 

The phase-field model appears to be a versatile method to study dendritic growth in 
general and sidebranching characteristics in particular. It has been shown to be adapted to 
take into account thermodynamical fluctuations p5|-p7|j, but alternative ways to introduce 
noise (such as the one employed in this work) appear to be appropriate to qualitatively 
reproduce the behaviour of sidebranching when the noise is of an external origin. 

Our results on sidebranching can also be helpful for further simulations where qualita- 
tively realistic sidebranching needs to be distinguished from that generated by numerical 
noise due to the round-off of the computer. 
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VII. FIGURE CAPTIONS 



FIG. 1. Amplitude vs. distance to the tip z of the sidebranching induced by internal noise 
in the experiments of Ref. [|18|]. Line: theoretical prediction of Eq. |2.10| with a overestimated 
value of S; crosses: experimental results (taken with permission from Ref. |18|). Quantities 
are expressed in micrometers. 

FIG. 2. Dendrite obtained with the phase-field model with a noise term included, as is 
described in the text. Ticks denote number of grid points. 

FIG. 3. Half-width h z (t) of the dendrite and its power spectrum P z (f) (P in square grid 
points times dimensionless time and frequency in cycles per dimensionless time) for z = 40 
grid points. 

FIG. 4. Half-width h z (t) of the dendrite and its power spectrum P z (f) (P in square grid 
points times dimensionless time and frequency in cycles per dimensionless time) for z = 100 
grid points. 

FIG. 5. Square root of the area under the spectral peak as a function of z (both in grid 
points). The line indicates exponential growth with z 2 ^ 5 . 
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